Stochastic Resonance in Two DimensionalLandau Ginzburg Equation 
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We study the mechanism of stochastic resonance in a two dimensional Landau Ginzburg equation 
perturbed by a white noise. We shortly review how to renormalize the equation in order to avoid 
ultraviolet divergences. Next we show that the renormalization amplifies the effect of the small 
periodic perturbation in the system. We finally argue that stochastic resonance can be used to 
highlight the effect of renormalization in spatially extended system with a bistable equilibria. 
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About twenty years asp, in studying the possible 
causes of the ice ages 00, we introduced the concept 
of stochastic resonance This process (hereafter, SR) 
requires that a physical system has its fundamental sym- 
metry broken in such a away that its phase space is not 
homogeneous. Thereby, the time spent in different re- 
gions of the phase space is not uniformly distributed, 
ff the system is stochastically forced, the probability of 
large fluctuation can be periodically locked to a small ex- 
ternal periodic forcing. The relevant implications of the 
SR are the generality of its set up and the limited band, 
in the parameter space, wherein the mechanism occurs. 
The latter allows for its exploiting in enhancing the sig- 
nal to noise ratio in a noisy system, while the former 
makes the concept applicables to a large class of sys- 
tems. Furthermore, we notice that the mechanism only 
requires that the system has large deviations properties 
and, therefore, the restriction to systems stochastically 
perturbed, is just a case amenable to explicit evaluation 
of the resonace range in the parameter space, fn practice, 
a stochastic forcing is required only for systems of gra- 
dient type (i.e., when a Lyapunov function exists surely) 
or whenever single point attractors exist. 

At the time of the publication above mentioned, we 
had doubts whether the mathematical process leading to 
stochastic resonance were a stochastic process at all. For- 
tunately in 4] has proved that a wide class of stochastic 
process obeying to a SR mechanism exists. Therefore we 
feel confident that, at least in a finite dimensional space, 
the SR is a well posed mathematical problem. When 
we consider stochastically perturbed infinite dimensional 
systems the problem may lead to a not consistent mathe- 
matics since the equations for the backward and the fore- 
ward probability density function may not exist. More- 
over, when a white noise in space and time is considered, 
the resulting equations are plagued by divergent corre- 
lation functions, fn a series of papers it was 
proved that these problems could be overcame by the 
method of the rcnormalisation theory. On the ground of 
these rigorous, mathematical results, in the present pa- 
per we wish to study the consequences of applying an 
external periodic (in time) forcing to a simple infinite 
dimensional stochastic equation. 



In particular, we discuss the SR phenomenon in two di- 
mensional Landau Ginzburg equation stochastically per- 
turbed by a white noise: 
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where 
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and is a white noise delta correletad both in space 
and in time. The field (j) is assumed to be real and defined 
on the torus [0, L]Ar[0, L],L = 27r, with periodic boundary 
conditions. Notice that the previous equation describes 
mean field theories nearby critical points, therefore, the 
use of white noise in space and time is mandatory when 
we do not want to introduce unknown physical scales. 
Moreover, unless m < 0, the theory has no broken sym- 
metries. 

As we mentioned, the phenomenon of SR has been 
widely studied in literature since it was introduced 
0,01 (for a review see ^3|). The study of SR in one space 
dimension, with a broken symmetry, has been presented 
m flU. fn this case, in fact, the stochastic equation does 
not show any particular novelty. As we shall see, the 
extension to more spatial dimension leads to interesting 
consequences. In particular we will consider their effects 
on the SR. 

As it has been discussed in that equation has 
no meaning in more than one space dimension unless it is 
renormalized. In D==2, the renomarlization amounts to 
change any non linear term in ^ by its Wick product, a 
rather well known procedure in quantum field theory. By 
using the Wick product, equations and (O becomes: 
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In eq. © we have introduced the cutoff field 4>a, where 
A is an ultraviolet cutoff in the Fourier space. Namely 
by denoting $(fc) the Fourier transform of (l){x) we have: 



0A 
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|fc|<A 
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and an analogous definition has been used for In 
© the quantity E„i{A) is the expectation value of sec- 
ond order moment of the field Z\ satisfying the linear 
stochastic differential equation: 



dZ\ = vAZj^ — mZty + e 

In two space dimensions one has: 

£2 vis? + 
i?™(A) = C-ln{ 



dt 



(5) 
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where C is a constant order one. Finally, let us recall 
that the quantity — 3Em{A)(j)\ is the Wick product 

One can show [TjI that eq. (|3Jl has a well defined limit 
for A going to infinity. Let us remark that this limit is 
achieved by keeping constant the finite volume. In this 
case one can show that equation ^ becomes: 

dt(j)A = v/\4>A — V e—jT- (7) 



dt 



where: 



(8) 



The value of p is determined by the following set of equa- 
tions: 



a 
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(9) 
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where (/)e is the equilibria of the effective potential, i.e. 
(j)e is determined by solving the equation: 



eff 







(11) 



The set of equations (|9I1UI11|) have a clear physical mean- 
ing shortly reviewed as follows. First of all, let us 
note that p is independent of A in the limit A — » cxd be- 
cause of eq. (jnil. Eq. is a redefinition of the Wick 
product in terms of a new constant a, which is defined in 
equation H10() as the curvature of the effective potential 
at one of its equilibrium values, namely </>£. The physical 
idea of eq. H1U|I is that, the counter term introduced in 
order to renormalize the Landau Ginzburg equation, pro- 
duces a double well potential. The position of the new 
minimum, determined by p, must be consistent with the 
fluctuations around the minimum itself. 

Using equations H9I10I11|I one gets a nonlinear equation 
for the variable x = namely 



X = -2 + 2Klogx 



(12) 



where K 



Let us remark, for future purpose, 
that it is possible to show that a solution of eq. H12|) 



always satisfies x > 2K. As discussed in the effective 
potential describes rather well the statistical properties 
of the space averge field of eq.(^. 

In general we can consider the variables $(^) defined 
as the average in space of 4>{x) on a box of side I. For 
I < L — 27r, the statistical properties of $(/) are no longer 
described by the same value of p because of the increas- 
ing effect, at small scales, of the stochastic fluctuations. 
In particular, for I — > 0, one finds that the statistical 
properties of $(^) are described by a single well potential 
centered around $(^) = 0. 

We want now to discuss the solutions of eq. ^ when a 
small periodic perturbation is added to the system and in 
particular we want to discuss the effect of renormalization 
of the mechanism on the stochastic resonance. 

To this aim, following the original discussion given in 
Q we consider the case of a time independent constant 
A added to the r.h.s. of The crucial point is the 
computation of the effective potential Vef / previously in- 
troduced. Equation (jSJ) are unchanged while equations 
(|10lll|l become: 



A = 



(13) 
(14) 



The effect of the constant perturbation is to change the 
equilibrium solution of the effective potential and there- 
fore to change the fluctuations around the equilibrium 
solution, which in turn change the value of p. The latter 
effect is due only to renormalization. 

In order to understand qualitatively the contribution 
due to renormalimation, one can compute a perturbative 
solution of the set of equations H9I13I14|I in power of A. 
At flrst order one obtains: 



Pi, 2 = PO + 
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(15) 
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where H — ^^J^^ ( = 2po) and all quantities with 
the index refer to values computed for ^ — 0. The 
previous definition of H implies H = < i because of 
the discussion on the solution of eq. Hl'ifl. Let us remark 
that we obtain two different values of p and (f>e because 
two different values of 4>eo are possible, one negative and 
the other one positive. Thus there exists two possible 
effective potential, each one describing the two well re- 
spectively. 

By using equations H15() and (|16() is now possible to 
compute the effective potential difference between the 
two stable solutions and the unstable solution. One gets: 
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In the limit of e — > we have iJ — > and the effect of 
renormalization disappears. For finite value of the noise 
and for g large enough the renormalization amplifies the 
asymmetry of the double well potential. Therefore, we 
should expect that, in two dimensions Landau Ginzburg 
equation, the mechanism of SR acts for smaller value 
of the forcing amplitude with respect to the one esti- 
mated by looking at the equilibrium probability distri- 
bution with no periodic forcing. 

In order to test the renormalization effect on SR, we 
have performed a numerical simulation on square lattice 
of 16x16 points. The parameters of our numerical sim- 
ulation are g = 220, v = 0.1, m = 0.1, e — 0.1. On the 
lattice, the quantity Emi^), defined in eq. lO, is given 
by: 




where 



2S{m,v,q,p) 



(19) FIG. 1: Power spectrum for the space average field ^{L) ob- 
tained for the numerical simulation of eq. Q. The peak 
corresponds to a period of 18 time units. 



^(m, V, q,p) — ma? + Av — 2vcos{2TTq/N) ~2vcos[2t:p / N) 

In equation (|19l wc introduce the lattice mesh a = 2ti/N . 
It follows that K — a^^. 

For different values of A, by applying Newton method, 
we have computed the solution of equations H9I13I14|I . 
using H19|l . It turns out that, with our numerical param- 
eters, po = 4.9. and Afi, ~ 2A, i.e, the external forcing is 
amplified by almost a factor 2. 

We have performed a long time integration of ^ with 
A = Q. The equilibrium probability distribution of ^{L) 
has been found to be bimodal with maxima centered at 
and the average transition time r between the 



and the average transition time r 

two equilibrium t = 9 in time units of m~^. Next, us- 
ing the analytical theory of SR discussed in Q , we have 
computed, by taking into account the amplification of 
the renormalization, the value of A for which the system 
should show a clear SR. We have found A = 0.1. Finally 
we have applied a periodic perturbation on the r.h.s. of 
eq. ^ with amphtude ^ = 0.1 and period equal to 
2t = 18, i.e. two times the average transition time. We 
have performed a long time integration corresponding to 
2048 time units. In figure we show the power spec- 
trum of the average field As one can clearly seen, 
there is an extremely well defined peak at frequency 113 
which corresponds to a period of 18 time units. 

In order to show the amplification effect on SR due to 
renormalization, we have numerically integrated the zero 
dimensional stochastic differential equation: 



dx ^ 

^ Pox- gx 



,2tn^ dW , , 



where A — Q.\ as before and the variance a of the white 
noise has been tuned in order to reproduce an average 
transition time of 9 (cr = 0.0165) for A = 0. With this 
choice of a we have been able to reproduce the equilib- 
rium probability distribuiont of $(L) at A = for the 



frequency 



FIG. 2: Same as in fig. lb for wero dimensional stochastic 
differential equation H2fl^ . Note that the amphtude of the 
peak at period of 18 time units is much smaller with respect 
to figure Q 



2D Landau Ginzburg equation. As for the previous case, 
we have integrated eq. (|20|l for 2048 time units. In figure 
we show the power spectrum of x. At variance with 
figure only a small effect of the periodic forcing is 
felt by the system, i.e. the system does not exhibits any 
SR. 

The comparison between figures and 10), indicates 
clearly that the renormalization effect is acting as an am- 
plifier of the external forcing, as previously discussed in 
the framework of a perturbation theory. 

It is quite interesting to compute the statistics of the 
transition times between the minima at different scales as 
described by the variables $(0- In absence of forcing the 
probability distribution of the transition time for '^(V) is 
exponential ,12] . In the stochastic resonance the proba- 
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time 

FIG. 3: Probability distribution of the transition time for 
$(!/), continues curve, and <E>(y|), continuos curve with sym- 
bols. 

bility distribution of the transition time for $(£) should 
show a well defined maximum at 9 time units. On the 
other hand, if we consider the transition time for $(/) 
with I <C L, then we expect that the forcing is not able 
to produce any SR. This is due to the fact that the pa- 



rameter describing the effective potential for small I are 
changed and, in particular, the average transition time 
between the minima becomes much smaller This 
effect is indeed observed in our numerical simulation. In 
fig. (2) we show the probability distribution of the tran- 
sition time for and ^(^)- As predicted, for the 
small scale, the probability distribution of the transition 
time is exponential. We can therefore reach the conclu- 
sion that, becuase of the renormalization, the SR mech- 
anism becomes scale depedent, a feature not previously 
observed. 

In this letter we have shown that renormalization can 
amplify the SR in two dimensional Landau Ginzburg 
equation. We want to remark that this effect can also be 
used to test, in a given experimental situation, whether or 
not a renormalization mechanism is acting in the system. 
In this case SR can be used as a tool to measure quanti- 
tatively the effect, if any, of renormalization in spatially 
exnteded systems. We remark that the present effect 
may be not restricted to system of gradient type and, in 
principle, could be detected in spatially extended deter- 
ministic systems. 
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